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A unified general approach is presented for construction of solutions of the char- 
acteristic initial value problems for various integrable hyperbolic reductions of Ein- 
■ stein's equations for space-times with two commuting isometries in General Relativ- 

ity and in some string theory induced gravity models. In all cases the associated 



(N 



linear systems of similar structures are used and their fundamental solutions ad- 
mit an alternative representation by two "scattering" matrices of a simple analytical 

> 

• structure on the spectral plane. The condition of equivalence of these representations 

o : 

, leads to the linear "integral evolution equations" whose scalar kernels and right hand 

sides are determined completely by the initial data for the fields specified on the two 

m ' 

initial characteristics. If the initial data for the fields are given, all field components 

O - 

Qi. of the corresponding solution can be expressed in quadratures in terms of a unique 

(3JT), solution of these quasi - Fredholm integral evolution equations. 



I. INTRODUCTION 

The Einstein's equations in General Relativity as well as in some string theory induced 
gravity models long ago were found to be integrable in a number of physically important 
cases. It is well known, that if the space-time admits an Abelian two-dimensional isometry 
group and all components of metric and matter fields and their potentials are dependent on 
some two of the four space-time coordinates only, the Einstein's equations for gravitational 
fields in vacuum as well as the Einstein - Maxwell field equations for electrovacuum fields 
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have been reduced to integrable equations. The similar reductions of the Einstein - Maxwell - 
Weyl equations for gravitational, electromagnetic and massless two-component spinor fields, 
the Einstein's equations for the ideal fluid with p = e, generalized Einstein's equations 
which arise as the string background equations and include electromagnetic, dilaton and 
axion fields, and some others are also integrable. 

The structures of all these integrable reductions of Einstein's equations can be described 
in a unified manner [1] . For this type of associated linear systems the space of local solutions 
of the field equations can be mapped by the space of monodromy data of a fundamental 
solution of associated linear system on the spectral plane. These space of monodromy 
data consists of a finite set of coordinate independent and arbitrarily chosen functions of 
spectral parameter. For example, any vacuum local solution is characterized by two and any 
electrovacuum solution - by four functions of spectral parameter. These functions can be 
calculated (in principle, at least) for any given solution of the field equations, and for any 
choice of these functions there exists a corresponding solution of the field equations. Thus, 
similarly to the well known inverse scattering transform, the direct and inverse problems 
of such "monodromy transform" arise and each possess a unique solution. To find the 
monodromy data for any given solution (the direct problem), one has to find for this solution 
the corresponding solution of the associated linear system with spectral parameter. Solution 
of the inverse problem, i.e. a construction of the solution for any given monodromy data 
functions, is equivalent to solution of some linear singular integral equations whose scalar 
kernels and right hand sides are expressed completely in terms of given monodromy data 
functions [2], [3]. 

The described above construction of the monodromy transform suggests immediately 
some methods for construction of infinite hierarchies of particular exact solutions of these 
field equations [4], [5] and provides a principle opportunity for solution of the initial or 
boundary value problems for the reduced Einstein's equations [6]. Farther development 
of this approach was suggested in [7], where for vacuum and electrovacuum fields some 
new, quasi - Fredholm form of linear integral equations substituting in this approach the 
nonlinear dynamical field equations was found. The structure of these integral equations 
is completely determined by the initial data for the fields on the mentioned just above two 
characteristics. Thus, this approach became well adapted for solution of the characteristic 
initial value problems, but its essential restriction was concerned with the conjecture of local 
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analyticity of solutions. A straightforward generalization of this approach to the solutions 
with non-analytical behaviour of fields and its application to solution of the characteristic 
initial value problem for plane gravitational or gravitational and electromagnetic waves with 
distinct wavefronts colliding on the Minkowski background was presented in [8]. 

In this paper, a generalization of the monodromy transform approach for other hyperbolic 
integrable reductions of Einstein's equations in General Relativity and string theory induced 
gravity models is presented. 

II. INTEGRABLE HYPERBOLIC REDUCTIONS OF EINSTEIN'S EQUATIONS 

Each of the mentioned above integrable reductions of Einstein's equations possess the 
space-time metric of the form 

ds 2 = g^dafdx" + g ab (dx a + oo^dx^) (dx b + u b u dx u ) (1) 

where fj,,u, ... — 1,2 and a, b, . . . = 3, 4 and the metric components g^ v , g a b, u a ^ are functions 
of two space-time coordinates x^ only. The electromagnetic, Weyl massless spinor, dilaton 
and axion fields, which can exist in these space-times, possess respectively the potentials 

Mx^), ip A {x»), <t>{x»), a t (xn (2) 

which are independent of the coordinates x a . Here the indices i,j, k, . . . = 1, 2, 3, 4, while the 
spinor indices A, B, ... — 1,2.) It seems useful to remined here that the metric functions u a M 
should be nonzero, if a nontrivial Weyl spinor field presents in the space-time (1). However, 
uj a ^ should vanish, if the spinor field is absent (ip^ = 0) and all the remaining fields together 
with their potentials, like, for example, the Kinnersly self-dual matrix potential for vacuum 
gravitational field, are independent of x a . 

In the "hyperbolic" case, all fields are functions of time and one of spacial coordinates. 
For this case it is convenient to use instead of x^ some pair of null coordinates (u,v), such 
that g^dx^dx" = f(u,v)dudv. Among these null coordinates, determined up to arbitrary 
transformations u = u(u) and v = v(v), we chose a pair of coordinates defined as 

follows. As is well known, the Einstein's field equations for metrics (1) and matter fields (2) 
imply that the function a{u, v) defined as det = « 2 should satisfy a simple equation 
d u d v a = 0. This allows to introduce another function f3{u,v) determined up to a real 
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additive constant by the relations: d u (5 = d u a and d v (5 = —d v a. Then, we define (£,r)) as 

£(u) = f3(u, v) + a(u, v), r)(v) = (3(u, v) — a(u, v). 

These coordinates are very convenient for investigation of integrability properties of the 
reduced Einstein's equations. 

III. EQUIVALENT MATRIX PROBLEMS 

The integrable hyperbolic reductions of Einstein equations mentioned above can be rep- 
resented by an equivalent complex N x N - matrix problems (N = 2 and iV = 3 for pure 
vacuum and electrovacuum cases respectively and the same is with presence of the Weyl 
spinor field, N = 4 for a generalized gravity model with electromagnetic, axion and dilaton 
fields) formulated for the four unknown N x iV-matrix functions [1] 

where w is a free complex ( "spectral" ) parameter. These should satisfy the following three 
groups of conditions. The first group is a linear system for ^ with case dependent structures 
of canonical forms of its coefficients 



U(e^) = ^^)U(o)-F- 1 (e^) 

(3) 



2i(w-0^* = U(£,»7)* 

2i(w-r,)dr,^ = V^,r,)^ 

where the transformation matrices T and Q depend on the field components and their 
derivatives, while the expressions for U( ), V( ), given below, are constant or, at most, 
depend on the spinor field current functions. The second group of conditions implies the 
existence for Eq. (3) of a Hermitian integral with case dependent constant matrix ft: 

wJh = WoH 

where the components of W also include the components of the metric (1) and the potentials 
of the matter fields (2). The third group of conditions are the normalization conditions 
imposed at some chosen "initial" point £ = £ j V — Vo 011 the values of ^ and W(£ , ??cb w ) — 
W„H: 

*(fo, Vo, w) = I, WoH = ^(w - Po)n + Go (5) 



dW 
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where /3 = (£ + ^o)/2 and G is also the case dependent constant matrix. 

In Eq. (3) the canonical forms U( ), V( ) of the matrices U and V (up to a permutation 
of diagonal elements) are 



U^ 2 = diagM) 



yN=2 
V (0) 



diag (i, 0) 



U^ 3 = diag(* + a(O,0,0) 



V iV=3 
V (0) 



diag (i, i, 0, 0) 



tj7V=4 
U (0) 

V^f 4 = diag 0,0) 



ft 



N=2 




^ 1 0^ 



diag (1 + 6(77), 0,0) 

where a(£), b(rj) are the arbitrarily chosen spinor field current functions [9]. The constant 
matrices ft should be chosen in the form: 

/ o 1 o\ 

1 
-1 
\ -1 J 

The chosen simplest normalization of the field components at the initial point (£0,770) leads 
to the following expressions for G (ct = (£ — ?7o)/2): 



N=3 
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n 



N=4 



G%= 2 = diag{4a 2 , 4}, G V = 3 = diag{4a 2 , 4, 1}, G V = 4 = diag{4« 2 , 4a 2 , 4, 4} 



IV. THE ANALYTICAL STRUCTURE OF * ON THE SPECTRAL PLANE 

The solution of the matrix problem (3) - (5) can be reduced to a decoupling pair of linear 
integral (quasi- Fredholm) equations with scalar kernels. For this, similarly to electrovacuum 
case [7], [8], we introduce two matrix functions Vf r + (£,u>) and &_(r),w) called there as "in- 
states" 

*+(S,w) = *(f,77o, w), *-(ri,w) = *(S ,ri,w). (6) 

Each of them is a fundamental solution of the linear system of ordinary differential equations 
supplied with the unit matrix as the initial conditions: 

2t(w - 0^* + = U(£, 770) • *+ ( 2i{w - 77)9„*_ = V(£ , V) • *- 

*+(&,«>) = I (*-{r] ,w)=I 
Hence, each of these "in-states" is determined completely by the values of the field compo- 
nents given on the characteristic £ = £ or V — Vo- 

The components of are holomorphic on the spectral plane outside the cut consisting 
of two non-overlapping segments of the real axis L + and L_, which join the points w = £ 
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with w = £ and w = r] with w = r/ respectively. Besides that, is holomorphic outside 
L + and is holomorphic outside L_ and near these cuts their local structures are 

L + : * = $ + {Z,r j ,w)®k+(w) + M + (Z,ri,w) 

(7) 

*+ = 0o+(£, w) <g> k+(iw) + M 0+ (£, u>) 

L_ : * = 0_(£,?7, w) (g) k_(w) + M_(£,?7, w) 

~ (8) 
*_ = t/> _(ri,w) ®k_(w) + M _(r7,w) 

Here 0±, 0o± are iV-component column- vectors and k±(u>) are iV-component row- vectors for 
the vacuum and electrovacuum cases (or with possible presence of the Weyl spinor field), i.e. 
when N = 2,3. However for the mentioned above string gravity model with electromagnetic, 
axion and dilaton fields (N = 4) the "columns" ip±, 0o± an d "rows" k-t take the values in 
the Grassman manifild (^^(C), i.e. they are 4x2- and 2 x 4- matrix functions respectively. 
For any (analytical) local solution of reduced Einstein's equations the N x iV-matrices M±, 
M ± and the vectors or matrices k±(iu) are single- valued (holomorphic) functions of w 
near the corresponding cuts L + or L_, whereas if)±, 0o± are the only branching solutions. 
The coordinate independent vectors (matrices) k ± (w) are related with the monodromy data 
which characterize every solution [3],[1]. 

Father, for all cases, for any local solutions analytical at the initial point (£o>?7o)> the 
branching character of 0±, 0o± is solution independent: 

L+ ■■ 0+ = A; 1 v+, 0o+ = A; 1 o+ 

L- : 0_ = AlV-, 0o- = AlVo- 



A+ = W^, A^= 



^ - 6 y w -r] 



where + , 0o+ an d 0-, 0o- are holomorphic near L + and L_ respectively; the functions A + 
and A_ are holomprphic functions of w outside L + and L_ respectively and X±(w = oo) = 1. 
However, for the local solutions with non-analytical behaviour at the initial point (£0,^0) as ) 
for example, at the point of collision of waves with distinct wavefronts [8], the character of 
branching of + , 0o+ and ip _ can be more complicate and solution-dependent. Besides 
that, in the presence of the spinor field, the character of branching of ^ also changes and 
we have to substitute everywhere 

A+ -> A+e +, A_ -> A_e , 0+ = - / -, a =- -. 

2 Jl+ w-( 2 J L _ w-( 

The defined above fragments of the algebraic structures of * and on the cuts L± will 

be used farther in our construction. 
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V. THE LINEAR INTEGRAL "EVOLUTION EQUATIONS" 

To construct the general solutions for the matrix problems (3) - (5) we express Sf' in two 
alternative forms in terms of "scattering" matrices x± : 

* = V, w) ■ w), * = V, w) ■ *_(?7, w). (9) 

Using the methods described in [3], [7], we can show that X- is holomorphic outside L + , 
while x+ is holomorphic outside L_ and they can be represented as the Cauchy integrals 
over these cuts. In the integrands, the jumps of the "scattering" matrices on L± can be 
expressed in terms of the functions defined in (7), (8). Then the condition of identity of 
the alternative representations (9) for \I> gives rise to the basic system of quasi - Fredholm 
linear integral equations: 

4>-(£,V, r -) - / S -(V,i~-, (+)<!>+(£, V,(+)d(+ = 0o-(»7,t_) 
Jl+ 

where r + ,(+ G L + and t_,(- e ^-j the functions </>+(r+), 0_(r_) and their initial values 
o +(r + ), 4> -(r-) are the jumps [t/>+] r+ , [V>-] r _ and [^o+]r+, [^o-]t_ respectively. The 
expressions for the scalar kernels are 

q (t ~ r \ ( m -(^C-) • 4>o+(£,t+)) . , x v 

^+(?> r +>C-) = -77 v , m_ = k_(wj • W + 

SJ } = (m + (,.C + ).^0^-)) > m+ = k+(ttl) . 

Thus, the coefficients of our generalized integral "evolution equations" are completely de- 
termined by the characteristic initial data for the fields. In terms of their solution, all of 
the field components can be determined in quadratures. In particular, the N x iV-matrix 
functions R±(£,?7) defined by the expansions x± — l + (l/w)H± + 0(l/w 2 ) can be calculated 

as 

1 f ~ 1 f ~ 
R+ = : / V- <E> m_ d(-, R = : / ip + ® m + d (+. 



™ J L _ m JL+ 

In terms of these matrix functions, for solution of (3) - (5) we obtain 

U(£, 77) = U(e, 770) + 2i^R+, V(£, 77) = V(&>, r/) + 2id„R_ 

and therefore, the derived linear integral evolution equations solve the characteristic initial 
value problems for considered integrable hyperbolic reductions of Einstein's equations. 



8 



Acknowledgments 

This work is supported partly by the Russian Foundation for Basic Research (grants 02- 
01-00729 and 02-02-17372) and by the Russian Academy of Sciences (the program "Nonlinear 
Dynamics and Solitons") 



[1] G. A. Alekseev, in Proc. of the Ws. " Nonlinearity, Integrability and all that: Twenty years 

after NEEDS'79", (World Scientific, Singapore, 2000); gr-qc/9911045. 

[2] G. A. Alekseev, Sov. Phys. Dokl. (USA) 30, 565 (1985). 

[3] G. A. Alekseev, Proc. of the Steklov Inst. Maths. 3, 215 (1988). 

[4] G. A. Alekseev and A. A. Garcia, Phys. Rev. D53, 1853 (1996). 

[5] G. A. Alekseev and J. B. Griffiths, Phys. Rev. Lett. 84, 5247 (2000). 

[6] G. A. Alekseev, in Proc. of the Ws. NEEDS'92, (World Scientific, Singapore, 1993). 

[7] G. A. Alekseev, Theor. Math. Phys., 129, 1466 (2001). 

[8] G. A. Alekseev and J. B. Griffiths, Phys. Rev. Lett. 87, 221101 (2001). 

[9] G. A. Alekseev, Sov. Phys. Dokl. (USA) 28, 133 (1983). 



